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Abstract. A proper understanding of boundary- value problems is essential in the at- 
tempt of developing a quantum theory of gravity and of the birth of the universe. The 
present paper reviews these topics in light of recent developments in spectral geometry, i.e. 
heat-kernel asymptotics for the Laplacian in the presence of Dirichlet, or Robin, or mixed 
boundary conditions; completely gauge-invariant boundary conditions in Euclidean quan- 
tum gravity; local vs. non-local boundary-value problems in one-loop Euclidean quantum 
theory via path integrals. 
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1. Introduction 



The aim of theoretical physics is to provide a clear conceptual framework for the wide 
variety of natural phenomena, so that not only are we able to make accurate predictions to 
be checked against observations, but the underlying mathematical structures of the world 
we live in can be thoroughly understood by the scientific community. What are therefore 
the key elements of a mathematical description of the physical world? Can we derive all 
basic equations of theoretical physics from a few symmetry principles? What do they tell 
us about the origin and evolution of the universe? Why is gravitation so peculiar with 
respect to all other fundamental interactions? 

The above questions have received careful consideration over the last decades, and have 
led, in particular, to several approaches [1] to a theory aiming at achieving a synthesis of 
quantum physics on the one hand, and general relativity on the other hand. This remains, 
possibly, the most important task of theoretical physics. The need for a quantum theory 
of gravity is already suggested from singularity theorems in classical cosmology [2] . Such 
theorems prove that the Einstein theory of general relativity leads to the occurrence of 
space-time singularities in a generic way. At first sight one might be tempted to conclude 
that a breakdown of all physical laws occurred in the past, or that general relativity 
is severely incomplete, being unable to predict what came out of a singularity. It has 
been therefore pointed out that all these pathological features result from the attempt 
of using the Einstein theory well beyond its limit of validity, i.e. at energy scales where 
the fundamental theory is definitely more involved. General relativity might be therefore 
viewed as a low-energy limit of a richer theory [3] , which achieves the synthesis of both the 
basic principles of modern physics and the fundamental interactions in the form presently 
known. 

Within the framework just outlined it remains however true that the various ap- 
proaches to quantum gravity developed so far suffer from mathematical inconsistencies, 
or incompleteness in their ability of accounting for some basic features of the laws of na- 
ture. From the point of view of general principles, the space-time approach to quantum 
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mechanics [4] and quantum field theory [5] , and its apphcation to the quantization of grav- 
itational interactions [6-9], remains indeed of fundamental importance. When one tries 
to implement the Feynman sum over histories one discovers that, already at the level of 
non-relativistic quantum mechanics, a well defined mathematical formulation is only ob- 
tained upon considering a heat-equation problem. The measure occurring in the Feynman 
representation of the Green kernel is then meaningful, and the propagation amplitude 
of quantum mechanics in flat Minkowski space-time is obtained by analytic continuation. 
This is a clear indication that quantum-mechanical problems via path integrals are well un- 
derstood only if the heat-equation counterpart is mathematically well posed. In quantum 
field theory one then deals with the Euclidean approach, and its application to quantum 
gravity [10] relies heavily on the theory of elliptic operators on Riemannian manifolds [11- 
14] . To obtain a complete picture one has then to specify the boundary conditions of the 
theory, i.e. the class of Riemannian geometries with their topologies involved in the sum, 
and the form of boundary data assigned on the bounding surfaces [15-29]. 

In particular, recent work [30,31] has shown that the only set of local boundary con- 
ditions on metric perturbations which are completely invariant under infinitesimal diffeo- 
morphisms is incompatible with a good elliptic theory. More precisely, while the resulting 
operator on metric perturbations can be made of Laplace type and elliptic in the interior 
of the Riemannian manifold under consideration, the property of strong ellipticity [12-14] 
is violated. This is a precise mathematical expression of the requirement that a unique 
smooth solution of the boundary-value problem should exist which vanishes at infinite 
geodesic distance from the boundary. This opens deep interpretive issues, since only for 
gravity does the requirement of complete gauge invariance of the boundary conditions turn 
out to be incompatible with a good elliptic theory. It is then impossible to make sense 
even just of the one-loop semiclassical approximation, because the functional trace of the 
heat operator is found to diverge (cf. section 6). 

We have been therefore led to consider, in our more recent research, non-local bound- 
ary conditions for the quantized gravitational field at one-loop level [32-35]. On the one 
hand, such a scheme already arises in simpler problems, i.e. the quantum theory of a free 
particle subject to non-local boundary data on a circle [36]. One then finds two families 
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of eigenfunctions of the Hamiltonian: surface states which decrease exponentially as one 
moves away from the boundary, and bulk states which remain instead smooth and non- 
vanishing. The generalization to an Abelian gauge theory such as Maxwell theory can fulfill 
non-locality, ellipticity and complete gauge invariance of boundary conditions providing 
one learns to work with pseudo-differential operators in one-loop quantum theory [37]. 
On the other hand, in the application to quantum gravity, since the boundary operator 
acquires new kernels responsible for the pseudo-differential nature of the boundary-value 
problem, one might hope to be able to recover a good elliptic theory under a wider variety 
of conditions [32-35]. 

Boundary field theory is therefore highly relevant for understanding quantum cosmol- 
ogy, quantum gravity and the foundations of quantized gauge theories, and it has deep 
roots in global analysis and spectral geometry, since, after imposing a supplementary con- 
dition, the gauge-field operator in the path integral for such theories [5,9] is either of 
Laplace type or non- minimal (on working with positive-definite metrics). In agreement 
with our pedagogical aims, we first review operators of Laplace type in section 2, the func- 
torial method for heat-kernel asymptotics in section 3 and the topic of mixed boundary 
conditions in section 4. Section 5 is then devoted to gauge-invariant boundary conditions 
for the quantized gravitational field, while recent developments and open problems are 
outlined in section 6. 

2. Asymptotics of the Laplacian on manifolds with boundary 

Following Branson and Gilkey [38], we are interested in a second-order differential 
operator, say P, with leading symbol given by the metric tensor (more precisely, with scalar 
leading symbol) on a compact m-dimensional Riemannian manifold M with boundary dM. 
Denoting by V the connection on the vector bundle V over M, our assumption implies 
that P, called an operator of Laplace type, reads 

P = -g^'^a^b - E, (2.1) 
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where E is an endomorphism of V (i.e. the 'potential' term in the physics-oriented htera- 
ture). The heat equation for the operator P is 



- + P)^ = 0. (2.2) 



By definition, the heat kernel is the solution, for t > 0, of the equation 



(^■^ + U{x,x';t) =0, xandx'eM, (2.3) 
subject to the boundary condition 

BU(x,x';t)] =0, (2.4) 

idM 

jointly with the (initial) condition 

lim / U{x,x';t)p{x')dx' = p{x), (2.5) 

which is a rigorous mathematical expression for the Dirac delta behaviour as t ^ 0+. The 
heat kernel can be written as 

U{x, x'; t) = J2 nn) i^)<P{n) {x')e-^'-^\ (2.6) 

(n) 

where {<^(„)(a;)} is a complete orthonormal set of eigenfunctions of P with eigenvalues 
A(„). The index n is enclosed in round brackets, to emphasize that, in general, a finite 
collection of integer labels occurs therein. 

Since, by construction, the heat kernel behaves as a distribution in the neighbour- 
hood of the boundary, it is convenient to introduce a smooth function, say / G C°°(M), 
and consider a slight generalization of the trace function (or integrated heat kernel), i.e. 

TrL2 (^fe~*^^ . It is precisely the consideration of / that makes it possible to deal properly 

with the distributional behaviour of the heat kernel near dM. A key idea is therefore to 
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work with arbitrary /, and then set / = 1 only when all coefficients in the asymptotic 
expansion 



Tr, 



oo 

~ (47rt)-^/2 J]r/2a^/2(/,P,S) (2.7) 

n=0 

have been evaluated. The term U {x, x; t) is called the heat-kernel diagonal. By virtue of 
Greiner's result [39], the heat- kernel coefficients a^/2(/, P, B), which are said to describe 
the global asymptotics, are obtained by integrating local formulae. More precisely, they 
admit a split into integrals over M (interior terms) and over dM (boundary terms). In 
such formulae, the integrands are linear combinations of all geometric invariants of the 
appropriate dimension (see below) which result from the Riemann curvature -R^^cd °f 
background, the extrinsic curvature of the boundary, the differential operator P (through 
the endomorphism E), and the boundary operator B (through the endomorphisms, or 
projection operators, or more general matrices occurring in it). With our notation, the 
indices a,b, ... range from 1 through m and index a local orthonormal frame for the tangent 
bundle of M, TM, while the indices range from 1 through m — 1 and index the 

orthonormal frame for the tangent bundle of the boundary, T{dM). The boundary is 
defined by the equations 

dM: ^ y^ix), (2.8) 

in terms of the functions y°'{x), being the coordinates on dM, and the those on M. 
Thus, the intrinsic metric, 7^ , on the boundary hypersurface dM, is given in terms of the 
metric gab on M by 

iij-9aby'',iy',j. (2.9) 

On inverting this equation one finds (here, n^" = N"- is the inward-pointing normal) 

g'^b ^ q'^b _^ n^n^ (2.10) 

where 

Q^' = y\j 1'' ■ (2.11) 



The tensor q"''' is equivalent to 7*-' and may be viewed as the induced metric on dM, in its 
contravariant form. The tensor q°'f^ is a projection operator, in that 

q\ <l\ = (2-12) 

n'' = 0. (2.13) 

The extrinsic-curvature tensor Kab (or second fundamental form of dM) is here defined 
by the projection of the covariant derivative of an extension of the inward, normal vector 
field n: 

Kab ^ n,,d q\ q\, (2.14) 

and is symmetric if the metric-compatible connection on M is torsion-free. Only its spatial 
components, Kij, are non- vanishing. 

The semicolon ; denotes multiple covariant differentiation with respect to the Levi- 
Civita connection Vm of M, while the stroke | denotes multiple covariant differentiation 
tangentially with respect to the Levi-Civita connection Vqm of the boundary. When 
sections of bundles built from V are involved, the semicolon means 

and the stroke means 

VaM® I+I® V. 

The curvature of the connection V on F is denoted by il. 

When Dirichlet or Robin boundary conditions are imposed on sections of V: 

[(f>]dM = 0, (2.15) 

or 

in^Va + S)(p\ =0, (2.16) 

L J dM 

the global asymptotics in (2.7) is expressed through some universal constants 

•{cKj, 6j, Cj, dj, Cj} 
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such that (here Rab = abc ^^i^^i tensor, R = and □ = VVa = g°-^VaS/b) 



ao{f,P,B)= f TV(/), 
Jm 

a,/2{f,P,B) = S{47ry/^ [ Trif), 

JdM 



a^{f,P,B) = - I Tr aifE + a2fR 



M 



JdM 



bof{trK) + bif,N + b2fS 



as/2{f,P,B) = -(47r)i/2 / Tr f {c^E + c^R + C2W 

yo JdM L 

+ f;N (c^itvK) + CgS^ 

a2{f,P,B) = ^J^Tr f (as\jE + a^RE + a^E^ 
+ ae\jR + c^rR^ + agRabR"'' 

+ ag RabcdR'^'^^'^ + Q^lO^ab^"^ j 



NiN 



+ 



Tr 



fidiE-N + d2R,N + t^3(trK)u ' + dAK, 



360 jQM 

+ d^E{irK) + deRitiK) + d^W j^n^iirK) + dsRiNjNK'^ 
+ dgR\ijK'^ + dio(trK)3 + diiKijK'^itiK) 

+ di5-Si2 + dieSR\^j^ + dnSitiKf + digSKijK'^ 
+ digS\tTK) + d2oS^ + d2iS^. I*) + f■,N[e^E + eai? 

+ esi^Viiv + e4(trK)2 + e^K^jK'^ + eg^ltriC) + eg^^) 
f-NN(eQ{iTK) + eio^j + 67/.^^ 



+ 



(2.17) 
(2.18) 



(2.19) 



(2.20) 



(2.21) 



These formulae may seem to be very complicated, but there is indeed a systematic way to 
write them down and then compute the universal constants. To begin, note that, if k is 
odd, aki2{f 1 P-,B) receives contributions from boundary terms only, whereas both interior 
terms and boundary terms contribute to aif/2{f,P,B), if k is even and positive. In the ai 
coefficient, the integrand in the interior term must be linear in the curvature, and hence it 
can only be a linear combination of the trace of the Ricci tensor, and of the endomorphism 
E in the differential operator. In the 02 coefficient, the integrand in the interior term 
must be quadratic in the curvature, and hence one needs a linear combination of the eight 
geometric invariants [6] 

□ Z? E? ZT" ZT'S I I p r?2 p Tjab p pabcd o O^^b 

In the ai coefficient, the integrand in the boundary term is a local expression given 
by a linear combination of all invariants linear (or of the same dimension as terms linear) 
in the extrinsic curvature: trKjS and f-N- In the 03/2 coefficient, the integrand in the 
boundary term must be quadratic in the extrinsic curvature. Thus, bearing in mind the 
Gauss-Codazzi equations, one finds the general result (2.20). Last, in the 02 coefficient, 
the integrand in the boundary term must be cubic in the extrinsic curvature. This leads 
to the boundary integral in (2.21), bearing in mind that f-N is on the same ground of a 
term linear in Kij, while f-NN is on the same ground of a term quadratic in Kij. Note 
that the interior invariants are built universally and polynomially from the metric tensor, 
its inverse, and the covariant derivatives of R^^^, ^ab and E. By virtue of Weyl's work on 
the invariants of the orthogonal group [12,38], these polynomials can be formed using only 
tensor products and contraction of tensor arguments. Here, the structure group is 0{m). 
However, when a boundary occurs, the boundary structure group is 0{m — 1). Weyl's 
theorem is used again to construct invariants as in the previous equations [38]. 

3. Functorial method 

Let T be a map which carries finite-dimensional vector spaces into finite-dimensional 
vector spaces. Thus, to every vector space V one has an associated vector space T{V). 
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The map T is said to be a continuous functor if, for all V and PF, the map 



T : Hom{V,W) — > Hom{T{V),T{W)) 

is continuous. 

In the theory of heat kernels, an application of the functorial method is the analysis 
of heat-equation asymptotics with respect to conformal variations. Indeed, the behaviour 
of classical and quantum field theories under conformal rescalings of the metric 

gab = ^'^ gab, (3.1) 

with fl a smooth function, is at the heart of many deep properties: light-cone structure, 
conformal curvature (i.e. the Weyl tensor), conformal-infinity techniques, massless free- 
field equations, twistor equation, twistor spaces, Hodge-star operator in four dimensions, 
conformal anomalies [26]. In the functorial method, one chooses Q, in the form 

n = e^f, (3.2) 

where £ is a real- valued parameter, and / e C°°{M) is the smooth function considered in 
section 2. One then deals with a one-parameter family of differential operators 

P(s) = e-^'^ P(0), (3.3) 

boundary operators 

B{e)^e-^f B{0), (3.4) 
connections on V, endomorphisms E{e) of V, and metrics 

9ab{e) = e'^'f gab{0). (3.5) 

For example, the form (3.2) of the conformal factor should be inserted into the general 
formulae which describe the transformation of Christoffel symbols under conformal rescal- 
ings: 

= + n-' (^5% + 5« - gtcg'^^d) ■ (3.6) 
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This makes it possible to obtain the conformal-variation formulae for the Riemann tensor 
R'Xcd. tensors involving the effect of Christoffel symbols. For the extrinsic- 

curvature tensor defined in Eq. (2.14) one finds 



Kab = ^Kab - n„Vbfi + Qab^in)^, 



(3.7) 



which implies 



(3.8) 



The application of these methods to heat-kernel asymptotics relies on the work by 
Branson and Gilkey [38]. Within this framework, a crucial role is played by the following 
'functorial' formulae {F being another smooth function): 



d 
de 



a,/2(l,e-2-^P(0) 



e=0 



(m-n)a„/2(/,P(0)), 



(3.9) 



d 
de 



a,/2(l,P(0)-£F) 



= a,/2_i(F,P(0)), 



e=0 



^a./2(e-^^/F,e-^^/p(0)) 



= 0. 



e=0 



(3.10) 



(3.11) 



Equation (3.11) is obtained when m = n + 2. These properties can be proved by (formal) 
differentiation, as follows. 

If the conformal variation of an operator of Laplace type reads 



P{£) = e-^'f P(0) - eF, 



(3.12) 



one finds 



d 
de 



Tr^^(e-*^(^)) 



Tr, 



e=0 



2t/P(0) + tp)e-*^(°) 



(3.13) 



Moreover, by virtue of the asymptotic expansion (2.7), one has (the numerical factors 
(477)""^/^ are omitted for simplicity, since they do not affect the form of Eqs. (3.9)-(3.11); 
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following Branson and Gilkey [38] one can, instead, absorb such factors into the definition 
of the coefficients a„/2(/, P)) 



d_ 
di 



T^L^ if e-'^'^A ~ -1 f;(m - t-i a^/^l/, P(0)). (3.14) 

\ / ^ n=0 



Thus, if F vanishes, Eqs. (3.13) and (3.14) lead to the result (3.9). By contrast, if / is set 
to zero, one has P{e) — -P(O) — sF^ which implies 

, -1 oo 

-I £=0 n=0 
oo 

= t— /2j]t'/2a,/2_i(F,P(0)), (3.15) 

1=2 

which leads in turn to Eq. (3.10). Last, to obtain the result (3.11), one considers the 
two-parameter conformal variation 

P(£, 7) = e-2^-^P(0) - ^e-^^^F. (3.16) 

Now in Eq. (3.9) we first replace n by 77. + 2, and then set m = n + 2. One then has, from 
Eq. (3.16): 

^an/2+i(l,P(£,7)) = 0. (3.17) 

Equation (3.17) can be differentiated with respect to 7, i.e. (see (3.10)) 

, , d d 



= 



^^^^an/2+i(l,P(e,7)) = ^^«n/2+i(l,e-^^^(P(0) -7P)) 



= ^an/2(e-'^^P,e-^^^P(0)), (3.18) 

and hence Eq. (3.11) is proved. 

To deal with Robin boundary conditions one needs another lemma, which is proved 
following again Branson and Gilkey. What we obtain is indeed a particular case of a 
more general property, which is proved in section 6 of Ref. [28]. Our starting point is 
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M, a compact, connected one-dimensional Riemannian manifold with boundary. In other 
words, one deals with the circle or with a closed interval. If 



b:C°°{M) — >?R: 



is a smooth, real-valued function, one can form the first-order operator 



A= — -b, 
ax 



(3.19) 



and its (formal) adjoint 



ax 



From these operators, one can form the second-order operators 



(3.20) 



dx"^ 



- 6x - b' 



(3.21) 



Do = AA^ = 



dx'^ 



+ b^-b' 



(3.22) 



where bx = For Di, Dirichlet boundary conditions are taken, while Robin boundary 



conditions are assumed for On defining 



Jx — , ) J XX — 7 9' 

ax ox^ 



one then finds the result 



(n-1) 



On/2(/,^l) - On/2(/, ^2) 



an/2-i[fxx + 2bfx,Diy (3.23) 



As a first step in the proof of (3.23), one takes a spectral resolution for Di, say {9,^, Aj^}, 
where 9,^ is the eigenfunction of Di belonging to the eigenvalue X,^: 



D\ 6u — Xn 9i,. 



(3.24) 
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Thus, differentiation with respect to t of the heat-kernel diagonal: 



(3.25) 



yields 



d 



-UiDux,x;t) = - ^A.e-*^" OUx) = - J] e"*^" e,)e,. (3.26) 



Moreover, for any 0, the set 



provides a spectral resolution of D2 on Ker(D2)"'"5 and one finds 

^^U{D2,x,x;t) = -Yl ^'(^) 



- e-'^^{Ae,){Ae,). 



(3.27) 



Bearing in mind that AO^, = if Aj^ = 0, summation may be performed over all values of 
z^, to find 

2^^\^{D^,x,x]t) -U{D2,x,x]t) 

= 2 J]e-*^'' (e'X - b'el - hHl) + (C - he^){e', - he,) 



e'X - b'el + {elf - 2beX 



(3.28) 



On the other hand, differentiation with respect to x yields 



- 2b^ U{D^ , X, x; t) = ^ e-*^-- (2eX - 2&^') , (3.29) 
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which imphes 



2/ 



d_ 

m 



U{Di,x,x;t)-U{D2,x,x;t)^. 



(3.30) 



We now integrate this formula over M and use the boundary conditions described before, 
jointly with integration by parts. All boundary terms are found to vanish, so that 



14 {r.-^»)u{nu.,.;t)d. 



-^ + 2b-^\U{Di,x,x;t)dx 
= 2/- [[/(A, X, x; t) - U{D2, X, x; t)] . 



(3.31) 



Bearing in mind the standard notation for heat-kernel traces, Eq. (3.31) may be re- 
expressed as 

„ d 



dt 



Tr 



/xx + 26/,')e-*^i 



(3.32) 



which leads to Eq. (3.23) by virtue of the asymptotic expansion (2.7). 

The algorithm resulting from Eq. (3.10) is sufficient to determine almost all interior 
terms in heat-kernel asymptotics. To appreciate this, notice that one is dealing with 
conformal variations which only affect the endomorphism of the operator P in (2.1). For 
example, on setting n = 2 in Eq. (3.10), one ends up by studying (the tilde symbol is now 
used for interior terms) 



ai(l,P)^^ / 

JM 



aiE + q:2-R 



= ai{E,R), 



(3.33) 
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which imphes 



ai(l, P(0) - eF) = &i{E, R) - ai{E - eF, R) 



= 1/ TV 

6 Jm 



ai{E -{E- eF)) + a2{R - R) 



Tr{aisF), 



M 



and hence 



d 



-ai(l,P(0)-eF) 



= J / IV(aiF) = ao(F,P(0))= / TV(F). 
5=0 " Jm Jm 



By comparison, Eq. (3.35) shows that 



CKl = 6. 



An analogous procedure leads to (see (2.21)) 

a2{l,P{0) - eF) = a2 P, Ric, Riem, ri) 

— a2 — eF, R, Ric, Riem, 



360 



Tr 



M 



a^ReF + a^{E'^ - {E - eF)'^) 



^ ( a^ReF + a^i - e^F^ + 2eEF) 
360 Jm I ^ r 



Now one can apply Eq. (3.10) when n = 4, to find 



^i^(l,P(0)-eF) 



= ai(F,P(0)) = ^ / Tr\aiFE + a2FR]. 

6 Jm L -I 

Equating the coefficients of the invariants occurring in the equation (3.38) one 



a2 = 



16 



as = 30ai = 180. 



(3.40) 



Furthermore, the consideration of Eq. (3.11) when n = 2 yields 



d 



-ai e 



de 

+ 2fTr{aiFE) 



r=0 6 Jm 



d 

de 



Tr{aiFE) 



e=0 



+ 



^Tr{a2FR) 
de 



+ 2fTiia2FR) } . 



e=0 



(3.41) 



At this stage, we need the conformal-variation formulae 



e=0 



-2fE+-im-2)\Jf, 



(3.42) 



d 
de 



R{e) 



-2fR-2{m-l)Uf- 



e=0 



Since we are studying the case m = n + 2 = 4, we find 



d_ 
de' 



aife-2^^F,e-2^^P(0))l ^\ j T,\{a^ - &a2)FUf 



0, 



(3.43) 



(3.44) 



which implies 



(y-2 = -^Oi\ = 1, 
6 



a4 = 60q;2 = 60. 



(3.45) 
(3.46) 



After considering the Laplacian acting on functions for a product manifold M = Mi x M2 
(this is another application of functorial methods), the complete set of coeflBcients for 
interior terms can be determined [38]: 



0:3 = 60 , CKe = 12 , CKy = 5 , ctg = —2 , ctg = 2 , aio = 30. 



(3.47) 



17 



As far as interior terms are concerned, one has to use Eqs. (3.9) and (3.11), jointly with 
two conformal-variation formulae which provide divergence terms that play an important 
role [38]: 

" d 



e=0 



A.,(^,e-/P(0)) 



(m-2)ai(/F,P(0)) 



(3.48) 



d 
de 



a2(P,e-2^^P(0)) 



-(m-4)a2(/P,P(0)) 



e=0 



6/' + 10 fR 



ab 



(3.49) 



Equations (3.48) and (3.49) are proved for manifolds without boundary in Lemma 4.2 
of Ref. [38]. They imply that, for manifolds with boundary, the right-hand side of Eq. 
(3.48), evaluated for P = 1, should be added to the left-hand side of Eq. (3.9) when n = 2. 
Similarly, the right-hand side of Eq. (3.49), evaluated for P = 1, should be added to the 
left-hand side of Eq. (3.9) when n = 4. Other useful formulae involving boundary effects 
are [38] 



Jm ^ ^ ;a JdM L ^ ^ 

U ' = f\i - (tri^)/;iv, 

/ □/ = -/ /;iV, 
JM JdM 



NiN 



NN 



(tri^)/;Ar]. 



(3.50) 
(3.51) 
(3.52) 

(3.53) 



IdM JdM 

On taking into account Eqs. (3.48)-(3.53), the application of Eq. (3.9) when n = 2,3,4: 
leads to 18 equations which are obtained by setting to zero the coefficients multiplying 

/;jv (when n = 2), 
18 



f;NN , /;jv(tri^) , f;NS (when n = 3), 

and 

f^^K^'^. , f\,n\ , f.,NS{tvK) , f.NS^ 

f;NNS , /|. S (when n = 4). 

The integrals of these 18 terms have a deep geometric nature in that they form a basis for 
the integral invariants. The resulting 18 equations are 

-bo{m - 1) - bi{m - 2) + ^62(m - 2) - (m - 4) = 0, (3.54) 

-co(m — 2) — 2ci(m — 1) + C2(m — 1) — C6(m — 3) = 0, (3.55) 
— -co{m — 2) + 2ci(m — 1) — C2 — 2cs{m — 1) 

-2c4-C5(m-3) + ;^C7(m-2) = 0, (3.56) 

-c-jim - 1) + C8(m - 2) - C9(m - 3) = 0, (3.57) 
-6(m - 6) + ^di(m - 2) - Idiim - 1) - e7(m - 4) = 0, (3.58) 

-60(m - 6) - 2di - d^{m - 1) + ^cii4(m - 2) - ei(m - 4) = 0, (3.59) 

z 

1 

-10(m - 6) - 2^2 - d^im - 1) + + T,^\birn - 2) - e2(m - 4) = 0, (3.60) 

z 

-drim - 1) - 4 + 2^9 - esim - 4) + ;^di6(m - 2) + 4(m - 6) = 0, (3.61) 

z 

l-d^im - 2) - 2d6{m - 1) + rf7(m - 1) + - e6(m - 4) = 0, (3.62) 
z 
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— -d^im — 2) + 2dQ{m — 1) — — dg — Sdio{m — 1) 

- 2c?ii + \dx'j{m - 2) - e4(m - 4) = 0, (3.63) 



— d% — d%{ra — — d\\{m — 1) — 3(ii2 

+ ^c?i8(m-2)-e5(m-4) = 0, (3.64) 

dzim - 4) - \d^{m - 2) + 2dQ{m - 1) - dr - dg - 4(m - 6) = 0, (3.65) 

d^im - 4) - rfg - dg{m - 3) + 4(m - 6) = 0, (3.66) 

(m - 4)di3 = 0, (3.67) 

— -zd\^{m — 2) + 2(ii5(m — 1) — d\% — 2dn{m — 1) 

Zt 

-2di8 + (ii9(m-2)-e8(m-4) = 0, (3.68) 

-dig{m - 1) + ^d2o(TO - 2) - e9(m - 4) = 0, (3.69) 

-diA{m - 2) - 2cii5(m - 1) + diQ{m - 1) - eio(m - 4) = 0, (3.70) 
2 

^di4(m - 2) - 2(ii5(m - 1) + diQ - d2i{m - 4) = 0. (3.71) 
z 

This set of algebraic equations among universal constants holds independently of the choice 
of Dirichlet or Robin boundary conditions. Another set of equations which hold for either 

Dirichlet or Robin boundary conditions is obtained by applying Eq. (3.11) when n = 

3,4. One then obtains five equations which result from setting to zero the coefficients 
multiplying 

^;iv/;7V (whcu U = 3) , 

f-NF-NN , f-NNF-N , f-NF-NitrK) , /-nF-nS (whenn = 4). 
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The explicit form of these equations is [38] 



-4C5 - 5C6 + -Cg = 0, 



-5e6 - 4e7 + 2eio = 0, 
2ei - 1062 + 563 - 267 = 0, 

— 26i + 1062 — 63 — IO64 — 265 — 566 + 667 + 268 = 0, 
— 568 + 469 ~ 56io = 0. 



(3.72) 

(3.73) 
(3.74) 
(3.75) 
(3.76) 



Last, one has to use the Lemma expressed by Eq. (3.23) when n = 2, 3, 4, bearing in 
mind that 

Ei = E{Dr) = -h^-h', (3.77) 



(3.78) 



For example, when n — 2, one finds 



ai (/, Di) - ai if, D2) - ao (^Ux + 26/^, L>i) = 0, 



(3.79) 



which implies (with Dirichlet conditions for Di and Robin conditions for D2) 



I [6/(^1 - E2) - - 125/J + / [ - b2fS - (3 + h)f,N 

Jm JdM 



0. 



(3.80) 



The integrand of the interior term in Eq. (3.80) may be re-expressed as a total divergence, 
and hence one gets 



/ \{l2-b2)fS+{3-h)f.N 

JdM L 



0, 



which leads to 



61 = 3 , &2 = 12. 



(3.81) 



(3.82) 



Further details concern only the repeated application of all these algorithms, and hence we 
refer the reader to Branson and Gilkey [38]. We should emphasize, however, that no proof 
exists, so far, that functorial methods lead to the complete calculation of all heat-kernel 
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coefficients. For the time being one can only say that, when Dirichlet or Robin boundary 
conditions, or a mixture of the two (see section 4) are imposed, functorial methods have 
been completely successful up to the evaluation of the 05/2 coefficient (see Ref. [27] and 
references therein). 

4. Mixed boundary conditions 

Mixed boundary conditions are found to occur naturally in the theory of fermionic 
fields, gauge fields and gravitation, in that some components of the field obey one set of 
boundary conditions, and the remaining part of the field obeys a complementary set of 
boundary conditions [25]. Here, we focus on some mathematical aspects of the problem. 
The framework of our investigation consists, as in section 3, of a compact Riemannian 
manifold, say M, with smooth boundary dM. Given a vector bundle V over M, we 
assume that V can be decomposed as the direct sum 

V = Vn(B Vd, (4.1) 

near dM. The corresponding projection operators are denoted by n„ and 11^, respectively. 
On Vn one takes Neumann boundary conditions modified by some endomorphism, say S, 
of Vn (see (2.16)), while Dirichlet boundary conditions hold on Vd- The (total) boundary 
operator reads therefore [12] 



Bf 

On defining 



© 




- dM 





(4.2) 

dM 



= n„ - Ud, (4.3) 



seven new universal constants are found to contribute to heat-kernel asymptotics for an 
operator of Laplace type, say P. In other words, the linear combination of projectors 
considered in Eq. (4.3) gives rise to seven new invariants in the calculation of heat-kernel 
coefficients up to a2'- one invariant contributes to a^/2{f,P,B), whereas the other six con- 
tribute to a2(/, P, B) (of course, the number of invariants is continuously increasing as one 
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considers higher-order heat-kernel coefficients) . The dependence on the boundary operator 
is emphasized by including it explicitly into the arguments of heat-kernel coefficients. One 
can thus write the general formulae (cf. (2.20) and (2.21)) 

a3/2(/,P,B) = ^(47r)i/2 / IyU/V'i, +a3/2(/,P,^), (4.4) 

»0 JdM L 



~a2if,P,B) 



1 



Tr 



360 



(4.5) 



where 03/2(7, P,B) is formally analogous to Eq. (2.20), but with some universal constants 
replaced by linear functions of '^^Uni^di and similarly for a2(/, P, B) and Eq. (2.21). 
The work in Refs. [12,40] has fixed the following values for the universal constants {/3i} 
occurring in Eqs. (4.4) and (4.5): 



A = -12, (32 = 60, ^3 = -12, P4 = -24, ^5 = -120, Pe = -18, Pr = 0. 



(4.6) 



To obtain this result, it is crucial to bear in mind that the correct functorial formula for 
the endomorphism S is 



d 
de 



S{e) 



£=0 



-fS+-{m-2)f.Mllr 



(4.7) 



This result was first obtained in Ref. [40], where the author pointed out that n„ should 
be included, since the variation of 5' should compensate another suitable variation only 
on the subspace Vn. The unfortunate omission of led to incorrect results in physical 
applications, which were later corrected by Moss and Poletti [41], hence confirming the 
analytic results in Refs. [18,21]. 
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5. Gauge-invariant boundary conditions for the gravitational field 

For gauge fields and gravitation, the boundary conditions are mixed, in that some 
components of the field (more precisely, a one-form or a two-form) obey a set of boundary 
conditions, and the remaining part of the field obeys another set of boundary conditions. 
Moreover, the boundary conditions are invariant under local gauge transformations pro- 
vided that suitable boundary conditions are imposed on the corresponding ghost zero-form 
or one-form. 

We are here interested in the derivation of mixed boundary conditions for Euclidean 
quantum gravity. The knowledge of the classical variational problem, and the principle 
of gauge invariance, are enough to lead to a highly non-trivial boundary-value problem. 
Indeed, it is by now well-known that, if one fixes the three-metric at the boundary in 
general relativity, the corresponding variational problem is well-posed and leads to the 
Einstein equations, providing the Einstein-Hilbert action is supplemented by a boundary 
term whose integrand is proportional to the trace of the second fundamental form [42] . In 
the corresponding 'quantum' boundary- value problem, which is relevant for the one- loop 
approximation in quantum gravity, the perturbations hij of the induced three-metric are 
set to zero at the boundary. Moreover, the whole set of metric perturbations hf^i, are 
subject to the so-called infinitesimal gauge transformations 



where V is the Levi-Civita connection of the background four-geometry with metric g, and 
(fii^dx'^ is the ghost one-form. In geometric language, the difference between h^i, and h^i, 
is given by the Lie derivative along (p of the four-metric g. 
For problems with boundary, Eq. (5.1) implies that 



where the stroke denotes, as usual, three-dimensional covariant differentiation tangentially 
with respect to the intrinsic Levi-Civita connection of the boundary, while Kij is the 
extrinsic-curvature tensor of the boundary. Of course, ipo and ipi are the normal and 



(5.1) 



hij = hij -h -I- KijifQ 



(5.2) 



24 



tangential components of the ghost one-form, respectively. Note that boundaries make 
it necessary to perform a 3+1 split of space-time geometry and physical fields. As such, 
they introduce non-covariant elements in the analysis of problems relevant for quantum 
gravity. This seems to be an unavoidable feature, although the boundary conditions may 
be written in tensor language. 

In the light of (5.2), the boundary conditions 



h 



= 



dM 



(5.3a) 



arc gauge invariant, i.e. 



h. 



dM 



0, 



(5.36) 



if and only if the whole ghost one-form obeys homogeneous Dirichlet conditions, so that 



^0 



J dM 



= 0, 



(5.4) 



= 0. 



dM 



(5.5) 



The conditions (5.4) and (5.5) are necessary and sufficient since </?o and ipi are independent, 
and three-dimensional covariant differentiation commutes with the operation of restriction 
at the boundary. Indeed, we are assuming that the boundary is smooth and not totally 
geodesic, i.e. Kij ^ 0. However, at those points of dM where the extrinsic-curvature 
tensor vanishes, the condition (5.4) is no longer necessary. 

The problem now arises to impose boundary conditions on the remaining set of metric 
perturbations. The key point is to make sure that the invariance of such boundary condi- 
tions under the infinitesimal transformations (5.1) is again guaranteed by (5.4) and (5.5), 
since otherwise one would obtain incompatible sets of boundary conditions on the ghost 
one-form. Indeed, on using the Faddeev-Popov formalism for the amplitudes of quantum 
gravity, it is necessary to use a gauge-averaging term in the Euclidean action, of the form 



327:Ga Jm 



(5.6) 
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where is any gauge-averaging functional which leads to self-adjoint elliptic operators 
on metric and ghost perturbations. One then finds that 

5^^{h) = ^^{h) - ^^{h) = .F/ (5.7) 

where .F^'^ is an elliptic operator that acts linearly on the ghost one-form. Thus, if one 
imposes the boundary conditions 

^^^ih)] =0, (5.8) 

L J dM 

and if one assumes that the ghost field can be expanded in a complete orthonormal set of 
eigenfunctions u[^^ of J^^^ which vanish at the boundary, i.e. 

V4'^ = H'^ (5-9) 

= (5.10) 
= 0, (5.11) 



dM 



the boundary conditions (5.8) are automatically gauge-invariant under the Dirichlet con- 
ditions (5.4) and (5.5) on the ghost. 

At a deeper level, the boundary conditions (5.3)-(5.5) and (5.8) are invariant under 

Becchi-Rouet-Stora-Tyutin transformations [43], but their consideration has been aban- 
doned after the proof of resulting lack of strong ellipticity in Refs. [30,31]. 

6. Recent developments and open problems 

The boundary conditions involving tangential derivatives of the field have been studied 
in detail not only in Refs. [25,28,30,31] but also in Refs. [44-46]. In particular, the 
work in Refs. [45,46] has produced valuable results in hcat-kcrnel asymptotics which are 
very important to get a geometric understanding of one-loop divergences in quantum field 
theory. 
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As far as Euclidean quantum gravity is concerned, we can see three main alternatives 
if it is approached from the point of view of spectral geometry (for a modern perspective 
on yet other issues the reader is referred to the paper by Vassilevich in this volume): 

(i) If one insists on using the completely gauge-invariant boundary conditions of section 
5 when a supplementary (i.e. gauge-fixing) condition of de Bonder type is used [47], one 
should prove that, since the trace of the heat semigroup can be then split into a regular 
part and a singular part [31], a regularized ^(0) value can be defined (cf. Ref. [48]) which 
is only aflfected by the regular part of Tri,2(e~*^). 

(ii) One can instead choose local and strongly elliptic boundary conditions for quantum 
gravity, which are, however, only partially gauge-invariant [49]. 

(iii) Last, but not least, one can resort to non-local field theory by considering a non- 
local gauge-fixing functional in the path integral. If one could prove that the resulting 
set of symbols for the gauge-field operator and boundary operator which are compatible 
with strong ellipticity of the boundary-value problem [32-35] is non-empty, one would 
understand from first principles why the universe starts in a quantum state but eventually 
reaches a classical regime, since non-local boundary conditions may lead to 'surface states' 
having precisely this property [36]. 

Moreover, it remains to be seen whether such investigations have an impact on current 
developments in string and brane theory [50]. 
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